Neutrino clustering in growing neutrino quintessence 
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A growing neutrino mass can stop the dynamical evolution of a dark energy scalar field, thus 
explaining the "why now" problem. We show that such models lead to a substantial neutrino 
clustering on the scales of superclusters. Nonlinear neutrino lumps form at redshift ?«1 and could 
partially drag the clustering of dark matter. If observed, large scale non-linear structures could be 
an indication for a new attractive "cosmon force" stronger than gravity. 
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Quintessence cosmologies could explain the observed 
order of magnitude of dark energy. This holds, in parti- 
cular, for cosmic attractors or scaling solutions, where the 
dark energy density decays in time as ~ t~ 2 . Such models 
generically predict the presence of an homogeneous dark 
energy fraction fL^ of similar size as the dark matter frac- 
tion fl m [H, EJ IE 0] ■ Nevertheless, upper bounds at high 
redshift restrict early dark energy to be £lj,(z > 5) < 0.1 
Q, while at present dark energy dominates f2^(£o) ~ 3/4. 
Realistic quintessence models need to explain the recent 
increase of (z) by a factor of around ten or more. It has 
been recently proposed [E 0] that a growing mass of the 
neutrinos may play a key role in stopping the dynamical 
evolution of the dark energy scalar field, the cosmon. For 
a slow evolution of the cosmon, the scalar potential acts 
like a cosmological constant, such that the equation of 
state of dark energy is close to w = — 1 and the expan- 
sion of the universe accelerates. In these models, the 
onset of accelerated expansion is triggered by neutrinos 
becoming non relativistic. For late cosmology, z ^ 5, the 
overall cosmology is very similar to the usual ACDM con- 
cordance model with a cosmological constant. 

An efficient stopping of the cosmon evolution by the 
relatively small energy density of neutrinos needs a 
cosmon-neutrino coupling that is somewhat larger than 
gravitational strength. This is similar to mass var ying 
neutrino models [| i, EE EI EI EE El EE EE EE EI 
"""E 22 1, even though the coupling in those models 
is generically much larger. In turn, the enhanced at- 



traction between neutrinos leads to an enhanced growth 
of neutrino fluctuations, once the neutrinos have become 
non-relativistic [E IE EE E3| ■ I n view of the small present 
neutrino mass, m v (to) < 2.3eV, and the time dependence 
of m Vl which makes the mass even smaller in the past, the 
time when neutrinos become non-relativistic is typically 
in the recent history of the universe, say zr w 5. Neu- 
trinos have been free streaming for z > zr, with a corre- 
spondingly large free streaming length. Fluctuations on 
length scales larger than the free streaming length are 
still present at zr, and they start growing for z < zr 
with a large growth rate. This opens the possibility that 
neutrinos form nonlinear lumps [E l2lj ] on supercluster 
scales, thus opening a window for observable effects of 
the growing neutrino scenario. 

In this letter we show that neutrino perturbations in- 



deed grow non- linear in these models. Non- linear neu- 
trino structures form at redshift z « 1 on the scale of 
superclusters and beyond. One may assume that these 
structures later turn into bound neutrino lumps of the 
type discussed in 2l|. The cosmon field within the neu- 
trino lumps does not vanish, and we see in our investiga- 
tion how the cosmon-fiuctuations are dragged by the neu- 
trino fluctuations. We compute the neutrino clustering as 
a function of scale and redshift, pointing out also how the 
growth of cold dark matter fluctuations is affected within 
these scenarios. Our investigation is limited, however, to 
linear perturbations. We can therefore provide a reliable 
estimate for the time when the first fluctuations become 
non-linear. For later times, it should only be used to give 
qualitative limits. 

A crucial ingredient in this model is the dependence of 
the neutrino mass on the cosmon field </>, as encoded in 
the dimensionless cosmon-neutrino coupling /?(</>), 
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with time 



and (3 < the neutrino mass increases 
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where m„ is a constant and (3 = (3 + dB/dlncj). The 
coupling P(<f>) can be either a constant [6| or, in gen- 
eral, a function of <f>, as proposed in [tJ within a particle 
physics model. The cosmon field <j) is normalized in units 
of the reduced Planck mass M= (87rG JV )- 1 / 2 ,and^~ 1 
corresponds to a cosmon mediated interaction for neutri- 
nos with gravitational strength. For a given cosmological 
model with a given time dependence of <j>, one can deter- 
mine the time dependence of the neutrino mass m v {t). 
For three degenerate neutrinos the present value of the 
neutrino mass m v {to) can be related to the energy frac- 
tion in neutrinos 
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The dynamics of the dark energy scalar field can be 
inferred from the Klein Gordon equation, now including 
an extra source due to the neutrino coupling, 
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with p v and p v — w v p v the energy density and pressure 
of the neutrinos. We choose an exponential potential 



V{4>) = M 2 U{4>) = M 4 e~ Q< ^ , 



(5) 



where the constant a is one of the free parameters of our 
model. 

The homogeneous energy density and pressure of the 
scalar field <j) are defined in the usual way as 

P* = ^2+ V i<t>) , V4, = £s-V(4>) , = & . (6) 
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Finally, we can express the conservation equations for 
dark energy a nd g rowing matter in the form of neutrinos 
as follows 22, 2311 : 



-3W(1- 
p' v = -3W(1 



3w v )p v 
- 3w v )p v 



(7) 



The sum of the energy momentum tensors for neutri- 
nos and the cosmon is conserved, but not the separate 
parts. We neglect a possible cosmon coupling to Cold 
Dark Matter (CDM), so that p' c = -ZHp c - 

For a given potential ([5]) the evolution equations for 
the different species can be numerically integrated, giv- 
ing the background evolution shown in FIGfT] (for con- 
stant pi) Q. The initial pattern is a typical early dark 
energy model, since neutrinos are still relativistic and 
almost massless. Radiation dominates until matter ra- 
diation equality, when CDM takes over. Dark energy is 
still subdominant and falls into the attractor provided 
by the exponential potential (see 22, [23[ for details). 
As the mass of the neutrinos increases with time, the 
term ~ /3p„ in the evolution equation for the cosmon (|4|) 
(or ([7|)) starts to play a more significant role, kicking <p 
out of the attractor as soon as neutrinos become non- 
relativistic. This resembles the effect of the coupled dark 
matter component in 2J|. Subsequently, small decaying 
oscillations characterize the <f> — v coupled fluid and the 
two components reach almost constant values. The va- 
lues of the energy densities today are in agreement with 
observations, once the precise crossing time for the end 
of the scaling solution has been fixed by an appropriate 
choice of the coupling p. At present the neutrinos are 
still subdominant with respect to CDM, though in the 
future they will take the lead (see @ for details on the 
future attractor solution for constant p). 

The evolution equations for linear perturbations (in 
Fourier space), in Newtonian gauge (in which the non 
diagonal metric perturbations are fixed to zero) pBj , read 
for the growing neutrino scenario: 

5' = 3W(w - c!)<^ 
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FIG. 1: Energy densities of neutrinos (dashed), cold dark 
matter (solid), dark energy (dotted) and photons (long 
dashed) are plotted vs redshift. For all plots we take a 
constant /3 = —52, with a — 10 and large neutrino mass 
m v = 2.11 eV. 
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8' v = 3(H - P{4>W){wu - c 2 v )8 v 
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The equations for the density contrasts 8i(k) = 
y J Si(x)exp(— ik ■ x)d 3 x (defined as the Fourier trans- 
formation of the local density perturbation <$j(x) = 
5 pi(x) / pi(x) over a volume V) involve the velocity per- 
turbations, which evolve according to 

v'x = -H(l - 3uuW - p((j))cf)'(l - 3w„)— v$ 

p$ 

5, 
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FIG. 2: Longitudinal density perturbation for CDM (solid), 
v (dashed) and (j> (dot-dashed) vs redshift for k — O.lh/Mpc 
(upper panel) and k — l.lh/Mpc (lower panel, A = 8Mpc). 
The neutrino equation of state (dotted) is also shown. The 
long dashed line is the reference ACDM. 
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As usual, the gravitational potentials obey 
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where ttto. is the anisotropic stress for the species a and 
the sound velocities are defined by cf = Spi/Spi. The 
perturbed pressure for <j> is 



$P<t> = -oH' o 



(14) 



a* a* 

and the anisotropic stress ttt^ = 0, as in uncoupled 
quintessence, since the coupling is treated as an external 
source in the Einstein equations. The linear perturbation 
of the cosmon, 5(f>, is related to via 

6(j) = fi'v^/k , 



-m<t>' 
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(15) 



+a 2 (3{(f>)(p ll - 3p„)] . 

Note that 5(f> can equivalcntly be obtained as the solution 
of the perturbed Klein Gordon equation: 
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The evolution of neutrinos requires solving the Boltz- 
mann equation in the case in which an interaction be- 
tween neutrinos and the cosmon is present [29( ■ The first 
order Boltzmann equation written in Newtonian gauge 
reads 1281 



^ + i2(k.n)* 

or e 

= i-(k • n)fc 
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(17) 



q 2 d(j) d In q 

where ty ps is defined as the perturbed term in the phase 
space distribution [H, [2!| : 

(18) 



f(x\T,q,nj) = f (q) [l + ^f ps (x\T,q,rij)] 



& and are the metric perturbations, x % the spatial co- 
ordinates, t is the conformal time, q = ap = qh is the 
comoving 3-momentum, e = e(4>) = ^/ q 2 + m„((/>) 2 a 2 , / 
is the phase space distribution and /o its zeroth-order 
term (Fermi-Dirac distribution). 

The Boltzmann hierarchy for neutrinos, obtained ex- 
panding the perturbation ^ ps in a Legendre series can 
be written in Newtonian gauge as 

... 1+ ^' rfln/o 

e 



Ks,0 = 



qk 



ps,l 



qk 
(2l + l)e 



d\nq 

x ek ^ d\n f Q 
dq dinq 



(19) 
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where [H 



1 q a 2 m 2 d In m v d In /o 
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This allows us to calculate the perturbed energy and pres- 
sure as well as the shear for neutrinos: 



Sp u = a 



q 2 fo(q) 



de 



OT 4 f 4 

Spy = — I -^fo(q) 



ps,0 



dcj) C 



dqdQ 
dqdQ , 
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(pu +p v )<Tv = -y a -4 / q 2 dq — /o(g)* ps ,2 ■ (21) 



The anisotropic stress is related to the shear via ttt v — 

3 
2p,. 



■^-(pv +Pv)cr and in our case 
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Note also that the unperturbed neutrino density and 
pressure read 



Pu = a 4 / q 2 dqdfle((t>)f (q) , 



p„ = -a 4 / q 2 dqdVl — {(j))} a {q) 



(23) 
(24) 



We numerically compute the linear density perturba- 
tions both using a modified version of CMBEASY 13 ll 
and, independently, a modified version of CAMB [34j |. 
written in synchronous gauge. The initial conditions are 
chosen such that S v = 6 7 = \f~A at redshift zi s ~ 1100, 
with A the primordial power spectrum amplitude as de- 
termined from the CMB anisotropies [131 . We plot the 
density fluctuations Si as a function of redshift for a fixed 
k in FIGOl The neutrino equation of state is also shown, 
starting from 1 /3 when neutrinos are relativistic and then 
decreasing to its present value when neutrinos become 
non relativistic. The turning point marks the time at 
which neutrino perturbations start to increase. At the 
scale of k = O.lh/Mpc (corresponding to superclusters 
scales) shown in FIGUJi, neutrino perturbations eventu- 
ally overtake CDM perturbations and even force (j> per- 
turbations to increase as well, in analogy with dark en- 
ergy clustering expected in 32, 33| within scalar tensor 
theories. Notice, however, that the scale at which neutri- 
nos form nonlinear clumps depends on the model param- 
eters, in particular the coupling (3, the potential param- 
eter a and the present days neutrino mass. Those are 
related to the neutrino free-streaming length, the range 
of the cosmon field and its mass. A detailed investigation 
of the parameter space will be performed in future work, 
but we mention that the model [3] with varying (3 gives 
qualitatively similar results. 
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FIG. 3: Longitudinal gravitational potential vs wavenumber 
k and values of redshift fixed to 0.5 (solid), 5 (dashed), 50 
(dot-dashed), in the linear approximation (note that the lines 
for the latter two redshifts overlap) . We also plot a reference 
ACDM model (dotted). 



We emphasize again that the linear approximation 
looses its quantitative reliability once one of the Si 
reaches one - for the wavelength shown in FIG [2] this 
first concerns neutrinos. Strong neutrino clumping could 
produce a gravitational potential that, in turn, drags the 
CDM fluctuations, as seen in the linear approximation 
in FIG1U However, once the neutrinos form strong non- 
linearities - neutrino lumps - one expects that nonlinear 
effects substantially slow down the increase of S v and even 
stop it. The magnitude of the CDM-dragging by neutri- 
nos is therefore not shown - it might be much smaller 
than visible in FIGHJ These remarks concern the quan- 
titative interpretation of all the following figures, which 
are always computed in the linear approximation. 

Nevertheless, the linear approximation demonstrates 
well the mechanisms at work. As can be seen from 
FIGO the gravitational potential shows a strong growth 
at the moment when neutrinos go nonlinear. Despite the 
small neutrino fraction strong neutrino fluctuations 
can become its main source. In turn, the gravitational 
potential will source the growth of CDM, whereas the 
cosmon field inhomogeneities are dragged by the even 
stronger interaction with the neutrinos. In the nonlinear 
regime, one expects that the neutrino structures decou- 
ple from the expansion and form stable lumps of the type 
described in [2l|. Typically, the gravitational potential 
in such static loops is much smaller than the huge val- 
ues obtained in the linear approximation in FIG(3J This 
demonstrates again that a quantitative understanding of 
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FIG. 4: Redshift of first non linearities vs the wavenumber 
k for CDM (solid), v (dashed) and cj> (dot-dashed). We also 
plot CDM for a reference ACDM model (dotted). 



FIG. 5: Longitudinal density perturbation for CDM and neu- 
trinos vs wavenumber k in the linear approximation, at red- 
shifts z = 0, 5. 



the gravitational potential and the CDM dragging has to 
wait for a better understanding of the nonlinear evolu- 
tion of neutrino fluctuations. This will also be crucial for 
an estimate of the role of neutrino clumping on the CMB 
anisotropics via the Integrated Sachs- Wolfe effect. 

In FIGH we plot the redshift z ni at which CDM, 
neutrinos and <fi become nonlinear as a function of the 
wavenumber k. The case of CDM in the concordance 
ACDM model with the same present value of £1^ and 
for massless neutrinos is also shown for reference (dot- 
ted line). The redshift z n i roughly measures when non- 
linearities first appear by evaluating the time at which 
8(z n i) = 1 for each species. The curves in FIGH] are ob- 
tained in the linear approximation, such that only the 
highest curves are quantitatively reliable. This concerns 
CDM for large k and neutrinos for small k. The sublead- 
ing components are influenced by dragging effects and 
may be, in reality, substantially lower. 

We can identify four regimes: i) At very big scales 
(larger than superclusters) the universe is homogeneous 
and perturbations are still linear today, ii) The range 
of length scales going from 14.5 Mpc to about 4.4 x 10 3 
Mpc appears to be highly affected by the neutrino cou- 
pling in growing matter scenarios: neutrino perturba- 
tions are the first ones to go nonlinear and neutrinos seem 
to form clumps in which then both the scalar field and 
CDM could fall into. Note that the effect of the neu- 
trino fluctuations on the gravitational potential induces 
CDM to cluster earlier with respect to the concordance 



ACDM model, where CDM is still linear at scales above 
~87 Mpc. iii) For lengths included in the range between 
0.9 Mpc and 14.5 Mpc, CDM takes over. That is in 
fact expected since neutrinos start to approach the free 
streaming scale. In this regime CDM drags neutrinos, 
and this effect may be overestimated in the linear ap- 
proximation. Notice that in our model CDM clusters 
later than it would do in ACDM. There are two rea- 
sons for this effect. At early times, the presence of a ho- 
mogeneous component of early dark energy, Q,^ ~ 3/a 2 , 
implies that Cl m is somewhat smaller than one and there- 
fore clustering is slower 35j. At later times, Q m is smaller 



than in the ACDM model since for the same fi^, part of 
1 — f2^ = Q m + VL V is now attributed to neutrinos. In con- 
sequence massive neutrinos reduce structure at smaller 
scales when they do not contribute to the clumping. The 
second effect is reduced for a smaller present day neutrino 
mass, iv) Finally, at very small scales (below clusters), 
CDM becomes highly non linear and neutrinos enter the 
free streaming regime, their perturbations do not growth 
and remain inside the linear regime. 

The longitudinal density perturbation for CDM and for 
neutrinos is shown in FIG(5] as a function of k, for two 
values of redshift z. For the present epoch (z = 0) the 
total size of the density contrast is presumably strongly 
overestimated in the linear approximation. Nevertheless, 
this figure visualizes the range of scales for which strong 
neutrino clumping is expected. In this range also the 
dragging effect on CDM is maximal. 
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Maximum neutrino clustering occurs on supercluster 
scales and one may ask about observable consequences. 
First of all, the neutrino clusters could have an imprint 
on the CMB-fluctuations. Taking the linear approxima- 
tion at face value, the ISW-effcct of the particular model 
presented here would be huge and strongly ruled out by 
observations. However, non-linear effects will substan- 
tially reduce the neutrino-generated gravitational poten- 
tial and the ISW-effcct. Further reduction is expected 
for smaller values of [3 (accompanied by smaller a). It 
is well conceivable that realistic models for the growing 
neutrino scenario lead to an ISW-effect in a range inter- 
esting for observations. Particular features are possible 
consequences of the oscillations in the neutrino sector. 
This typically leads to an ISW effect which can show 
structures as a function of the angular momentum I. 

A second possibility concerns the detection of nonlin- 
ear structures at very large length scales. Such structures 
can be found via their gravitational potential, indepen- 
dently of the question if neutrinos or CDM source the 
gravitational field. Very large nonlinear structures are 
extremely unlikely in the ACDM concordance model. 
An establishment of a population of such structures, 
and their possible direct correlation with the CMB-map 
a& I2E could therefore give a clear hint for 
"cosmological actors" beyond the ACDM model. For 
any flat primordial spectrum the gravitational force will 
be insufficient to produce large scale clumping, whatever 
the ingredients are. Large scale clumping could thus be 
an indication for a new attractive force stronger than 
gravity - in our model mediated by the cosmon. 
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